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We study the space-time correlation and response functions in nonequilibrium growth processes 
described by linear stochastic Langevin equations. Exploiting exclusively the existence of space 
and time dependent symmetries of the noiseless part of these equations, we derive expressions 
for the universal scaling functions of two-time quantities which are found to agree with the exact 
expressions obtained from the stochastic equations of motion. The usefulness of the space-time 
functions is illustrated through the investigation of two atomistic growth models, the Family model 
and the restricted Family model, which are shown to belong to a unique universality class in 1+1 
and in 2+1 space dimensions. This corrects earlier studies which claimed that in 2+1 dimensions 
the two models belong to different universality classes. 

PACS numbers: 81.15.Aa,05.70.Np,05.10.Gg 

I. INTRODUCTION 

Fluctuations arc omnipresent when analyzing surfaces and interfaces. These fluctuations can be equilibrium fluc- 
tuations, as encountered for example when looking at steps on surfaces, or they can be of nonequilibrium origin as 
it is the case in various growth processes. Well known examples of nonequilibrium surface fluctuations are found 
in kinetic roughening or nonequilibrium growth processes |l|, |2j, |3( as for example in thin film growth due to vapor 
deposition. Interestingly, both equilibrium and nonequilibrium interface and surface fluctuations can be described 
on a mesoscopic level through rather simple Langevin equations Q. In this approach the fast degrees of freedom 
are modeled by a noise term, thus yielding stochastic equations of motion for the slow degrees of freedom. In many 
instances the physics of dynamical processes is to a large extend captured by linearized Langevin equations [j| where 
one distinguishes whether the dynamics is purely diffusive or whether mass conservation has to be implemented. 

For purely diffusive dynamics (called model A dynamics in critical dynamics the linear Langevin equation 

can be written in the following way: 

^|^ = ^ 2 V 2 Mx,t) + A + ?? (x,t) (1) 

where /i(x, t) is the value of the macroscopic field h at site x at time t. In the physical context of fluctuating interfaces 
and growth processes in d+1 spatial dimensions h is the height field whereas x is the lateral position in the underlying 
d dimensional substrat lattice. In addition, vi > is the diffusion constant whereas A is the mean growth velocity 
(which may of course be zero). Finally, the random variable r/ models the noise due to the fast degrees of freedom. 
Depending on the physical problem at hand, either Gaussian white noise or spatially and/or temporally correlated 
noise is usually considered y|- 

In the context of kinetic roughening and nonequilibrium growth processes Eq. Q is called the Edwards- Wilkinson 
(EW) equation Q . This equation has been used for the description of many dynamical proces ses, as for example 
equilibrium step fluctuations with random attachment /detachment events at the step edge jlfl llll Il2j . This equation 
also describes the dynamics of a growing surface with a normal incidence of the incoming particles. An obliquely 
incident particle beam, however, generates anisotropics which can only be described by a more complex non-linear 
Langevin equation |l3j| . 

In growth processes with mass conservation the following linear Langevin equation (with v± > 0) 

<%(x,i) _ 4 



Of 



= -^ 4 V 4 /i(x,i) + A + ??(x,t) (2) 



has been proposed [Tl . Il5l ] . This equation is sometimes called the noisy Mullins-Hcrring (MH) equation. The noise 
term again reflects the physics of the investigated system. In the case of equilibrium fluctuations conserved noise must 
be considered, leading to the so-called model B dynamics 0,0- On the other hand, when studying out-of-equilibrium 
processes one can again focus on Gaussian white noise or on noise which is correlated in space and/or time [l6[ . 
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The Langevin equation (|2J) is used for example to describe film growth via molecular beam epitaxy 0, 0, , 
equilibrium fluctuations limited by step edge diffusion 0, 0] or even tumor growth [2(3, E| ■ 

An important notion in noncquilibrium growth processes is that of dynamical scaling. Dynamical scaling is nicely 
illustrated through the behavior of the mean-square width of the surface or interface which for a substrate of linear 
size L scales as |23 

W 2 (L,t) =L 2 <F{t/L z ) (3) 

where £ is the roughness exponent and z is the dynamical exponent. For EW we have z = 2 whereas for MH z = 4. 
The value of £ depends on whether correlated or uncorrelated noise is considered. 

Langevin equations have the drawback that they do not really mirror the atomistic processes underlying the 
fluctuations of the interfaces and surfaces. In order to capture the physics on the microscopic level one commonly 
designs simple atomistic models (characterized by some specific deposition and/or diffusion rules) which are then often 
studied numerically. However, it is not always clear what the corresponding Langevin equation is. Usually, numerical 
simulations are used in order to extract the exponents £ and z (see Eq. 0) which generally permit to relate the 
microscopic model to one of the Langevin equations (universality classes). These exponents, however, encode only 
partly the information given by a scaling behavior, as the scaling functions, like F(y) in Eq. 0, arc themselves 
different for different universality classes. 

In this paper we focus on two-point quantities as for example the space-time response and the space-time correlation 
functions which also display a dynamical scaling behavior. On a more fundamental level we show, using arguments first 
given in p3j , that in systems described by the equations and @ the scaling functions of these two-point quantities 
can be derived by exclusively exploiting the symmetry properties of the underlying noiseless, i.e. deterministic, 
equations. This approach, which is based on generalized, space and time dependent, symmetries of the dynamical 
system 123. Ep. has in the past already been applied successfully in the special case z = 2 to systems undergoing phase 
ordering |23tl26t I27I |28| and to nonequilibrium phase transitions j2{j . Here we show that local space-time symmetries 
also permit to fix (up to some numerical factors) the scaling functions of space-time response and correlation functions 
in cases where z = 4. On a more practical level we demonstrate the usefulness of the scaling functions of these two- 
point quantities (which depend on two different space-time points (x,i) and (y,s)) in the characterization of the 
universality classes of noncquilibrium growth processes. Whereas in the study of critical systems scaling functions, 
which are universal and characterize the different universality classes, are routinely investigated (and this both at 
equilibrium |30j and far from equilibrium |3ltl32| ). in nonequilibrium growth processes the focus usually lies on simple 
quantities like for example the exponents £ and z. There are some notable exceptions where scaling functions have 
been discussed (see, for example, |aE3): but these studies were in general restricted to one-time quantities. However, 
also in nonequilibrium growth processes scaling functions of two-point functions are universal and should therefore 
be very valuable in the determination of the universality class of a given microscopic model or experimental system. 
We illustrate this by computing through Monte Carlo simulations the two-point space-time correlation function for 
two microscopic models which have been proposed to belong to the same universality class as the Edwards- Wilkinson 
equation with Gaussian white noise [33, l34l l35l 13a . l37j . 

The paper is organized in the following way. In the next Section we discuss the EW and the MH equations in more 
detail and introduce the two-point functions. Section III is devoted to the computation of the exact expressions of the 
space-time correlation and response functions by Fourier transformation. These exact results show inter alia that the 
response of the system to the noise does not depend explicitely on the specific choice of the noise itself. In Section IV 
we discuss the space-time symmetries of the noiseless equations, whereas in Section V we show how these symmetries 
can be used for the derivation of the scaling functions of two-point functions. In Section VI we numerically study two 
microscopic growth models which have been proposed to belong to the Edwards- Wilkinson universality class. There 
has been a recent debate on the universality class of these models which we resolve by studying the scaling function 
of the space-time correlation function. Finally, in Section VII we give our conclusions. Some technical points are 
deferred to the Appendices. 



II. NOISE MODELIZATION AND SPACE-TIME QUANTITIES 



Our main interest in this paper is the investigation of space-time quantities in systems described by the quite general 
linear stochastic equations Q and . Setting the mean growth velocity A to zero (which can always be achieved by 
transforming into the co-moving frame) both cases can be captured by the single equation 

d -^ = -M-^) l h(K,t) + ^t) (4) 
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with / = 1 (EW) or I = 2 (MH). As it is well known, these equations of motion can be derived from a free field theory 


Depending on the physical context, different types of noise may be considered. For the EW case we shall discuss 
both Gaussian white noise (EW1) 

(r?(x, t)) = 0, (7 ? (x, t)r,(y, s)) = 2D5 d (x - y)6(t - s) (5) 

and spatially correlated noise (EW2) 

(77(x,t))=0, (r,(x,tMy,s))=2D\x-y\ 2p - d 5(t-s) (6) 

with < p < d/2. In the past, these two types of noise have been used in the modeling of noncquilibrium growth 
processes 0,0. If, however, one wishes to model thermal equilibrium interface fluctuations, as for example step 
fluctuations rate-limited by evaporation-condensation, one has to consider white noise with the Einstein relation 
D = i/jfcflT where T is the temperature and Ub the Boltzmann constant. For the MH case we also consider the 
noises © and ©, called MH1 and MH2 in the following. In this case, however, white noise can only be used in 
noncquilibrium situations as it breaks the conservation of mass encoded in the Langevin equation @. We shall not 
consider here the noisy Mullins- Herring equation with conserved noise which assures the relaxation towards equilibrium 
of a system with conserved dynamics, as this is covered by a forthcoming publication of one of the authors |3*sj . 

Two-time quantities have been shown in many circumstances to yield useful insights into the dynamical behavior 
of systems far from equilibrium [32| . Of special interest are space and time dependent functions as for example the 
space-time response i?(x, y, t, s) or the height-height space-time correlation 

C(x,y,t,s) = (h(x,t)h(y,s)) (7) 

where the brackets indicate an average over the realization of the noise. The space-time response, defined by 

S(h(x,t)) 



■R(x,y,*,s) 



(8) 



measures the response of the interface at time t and position x to a small perturbation j(y,s) at an earlier time s 
and at a different position y |3j|. For reasons of causality we have t > s. At the level of the Langevin equation the 
perturbation enters through the addition of j to the right hand side. Assuming spatial translation invariance in the 
directions parallel to the interface, we have 

C(x,y,t,s)=C(x.-y,t,s), R(x,y,t,s) = R(x- y,t,s). (9) 

The autocorrelation and autorcsponse functions are then defined by 

C(t,s) :=C(0,t,s), R(t,s) :=R(0,t,s). (10) 

It is well known that the systems discussed here present a simple dynamical scaling bevahiour (sec Eq. 101). For 
example, for the autorcsponse and the autocorrelation functions we have 

R(t, s) ~ s-^fnit/s), C(t, s) ~ s- b f c (t/s), (11) 

which defines the noncquilibrium exponents a and b. This terminology is well known from magnetic systems. Com- 
bining equations Q and Ullfl . one ends up with the scaling relation b = —2(/z, which relates b to the known 
exponents C and z. In addition the scaling functions Jr and fc define two additional exponents A# and Ac by their 
asymptotic behavior 

My) y ~°° y~ x " /z , fc(y) v ~°°y- Xc/z (12) 

where z is again the dynamical exponent introduced in Section I. Similarly, one obtains for the space-time quantities 
the following scaling forms: 

i?(x - y, t, s) ~ a—^flx - y\ z /a, t/s), C(x - y, t, s) ~ S - b F c (|x - y\ z /s, t/s). (13) 



III. RESPONSE AND CORRELATION FUNCTIONS: EXACT RESULTS 



This Section is devoted to the computation of space-time quantities by directly solving the Langevin equations 
and 121 in the physically relevant cases d = 1 and d = 2. These exact results will be used in the following in 
two different ways. In Section V we use these expressions in order to check whether our approach, which exploits 
exclusively the generalized space-time symmetries of the deterministic part of the equation of motion, yields the 
correct results. In addition, in Section VI we compare these expressions with the numerically determined scaling 
functions obtained for two different atomistic models in order to decide on the universality class of these models. 
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A. 2 = 2: The Edwards- Wilkinson case 



In order to compute the response of the surface/interface to a small perturbation we add the term j(x, t) to the 
right-hand side of the Langevin equation and then go to reciprocal space. In the EW case the solution of the resulting 
equation is (with d = 1, 2) 



h(k,t) 



-v 2 k 2 t 



dt , e»> kH '(f)(k,t')+j(Kt')) 



(14) 



where we denote by h(k, t), fj(k, t) resp. j(k, t) the Fourier transform of h(x,t), r](x,t) resp. j( x > t). We thereby 
assume flat initial conditions at time t = 0, h(x, 0) = 0, i.e. we prepare the system in an out-of-equilibrium state j4jj. 
This preparation enables us to study the approach to equilibrium for the EW1 case with a valid Einstein relation. 
For the corresponding study of equilibrium dynamical properties (as encountered in the recent experiments on step 
fluctuations [IollllLll2T| ) we have to prepare the system at t = — oo and replace in l|14|l the lower integration boundary 
by -co. We shall in the following concentrate on the out-of-equilibrium situation. 

Taking the functional derivative of (|14f) and transforming back to real space yields the result 



i?(x-y,M) 



dk 



gik(x-y) e -^2k (t-s) 



= r (t - s) 2 cxp - 



(x-y) 



4^2 (t 



(15) 
(16) 



with tq = p— /j=ja and t > s. The exponents a and \r as well as the scaling function can readily be obtained 
from the expression of the autoresponse function (see Equations l|ll|) and (|12|l 'l 



R(t,s) = r (<-s)~ 



(17) 



yielding a = X R = d and f R (y) ~ (y - l)a 1 . 

The expression for the space-time response is completely independent from choice of the noise term as long as 
(77(k, t)) = 0. This also holds for the MH case, as discussed in the next subsection. In Section V we shall discuss an 
alternative way of looking at this fact. 

A similar straightforward calculation yields for the space-time correlation the expression 



C(x-y,t,s) = / dt' 



dt" 



dk 

(2tt)° 



D ik-(x-y) -v 2 k 2 (t+s-t'-t") 



(f,(k,t')f,(-k,t")} 



(18) 



where we have exploited the spatial translation invariance of the noise correlator. For Gaussian white noise we then 
obtain for the space-time correlation 



C(x-y,i,s) =c |x-y 



1 



' 4z/ 2 (i + s) 



r 



i. 



y) 2 



4i/ 2 (t-s) 



(19) 



with Co = 2d+3 ^ d / 2 ■ The autocorrelation function for d ^ 2 [44| is obtained by using the known series expansion of 
the incomplete Gamma-functions T: 



C(t,s) 



2c (4^ 2 ) 1 - 
2-d 



s 



t 

- - 1 

s 



(20) 



from which we find b = | — 1, Ac = d and fc{y) ~ {y + I) 1 2 — {y — I) 1 5 • For the special case d = 2 one has to 
take the logarithmic behavior of the Gamma-functions into account which yields 



C(i, s) = cq In 



(21) 



For spatially correlated noise the space-time correlator can only be written as a series expansion: 



C(x - y, t, s) = £ (-1)" oW(p)|x - y| 2n \(t + a)-^ 



■2p+d-2)/2 



(*■ 



-(2n-2p+d-2)/2 



(22) 
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with a>a\p) = b W - 
one gets 



2 2p D r(p)r(ra-p+d/2) 



(2p-2n-d+2)Tr 2 - 3d / 2 r(d/2-p)u'2 



and b 



(1) 



3- ^ 



(2) 



(2n)> 



((2n)!!) 2 ' 



whereas for the autocorrclator 



C(t,s) = a<i?\p)s 1 -i + r 



t V-i+o ft 



(23) 



yielding b = f - 1 + p, X c = d - 2p and / c (y) ~ (y + ljM+P - (y - lf-t+P . 

Looking at the expressions l|19[) and l|22|) . we see that in both cases the space-time correlation has the following 
scaling form: 



C(x-y,M) = |x-y| Q F 



(x-y) 2 (x-y) 5 



(24) 



with a = 2 — d resp. 2 — d + 2p for EW1 rcsp. EW2. The scaling function F is then a function of the two scaling 
variables 



t+s 



and t J^ ■ In the non-equilibrium situation we discuss here the space-time correlation function is 



therefore not time translation invariant. For equilibrium systems which are prepared at t = —oo time translation 
invariance is of course recovered. It is worth noting that the scaling form l|24|l corrects the scaling forms given in Q 

where only a dependence on . was predicted far from equilibrium. 



B. z = 4: The Mullins-Herring case 



For the Mullins-Herring case we proceed along the same line as for the Edwards- Wilkinson case. We here only give 
the results for one-dimensional interfaces and refer the reader to the Appendix A for the two-dimensional case. The 
solution of the Langevin equation in Fourier space reads in the MH case 



dt'e 



(tKMO+jXM')) 



(25) 



o 



which after some algebra yields for one-dimensional interfaces the expression 

1 
"1 

y 4 



-1/4 



i?(x-y,M) = — i^Wi-s)) 

(x-y) 4 



5\ (1 3. (x-y) 4 
1 1 4j ^ 2 U'4'256^-s) 



256(^ 4 (t - s)) 



1/2 F (I 3 ( x -y) 4 

2 U' 2' 256^ 4 (t-s; 



(26) 



for the space-time response. Here the qF 2 functions are generalized hypergeometric functions. It is worth noting that 
exponentially growing contributions to the functions 0-F2 just cancel each other, yielding a response function which 
decreases for (x — y) 4 / (t — s) — ► 00, as it should. 

The autoresponsc function is straightforwardly found (with t > s): 



R(t,s) 



r(f) 



1/4 



(t- s y 



(27) 



which gives us the quantities a 



1 and f R (y) ~ (y - 1) 



space dimension d one has the relations a = 4 — 1, A# = c? and fn{y) 



One straightforwardly verifies that in any 

(y-i)" J 



As for the EW case we remark that the exact expressions ()26(l and Ij27(l are independent of the noise: we obtain the 
same results for Gaussian white noise and for spatially correlated noise. Let us add that we have here only considered 
perturbations which are not mass conserving. Whereas this is physically sound for the cases wc have in mind here, 
one usually considers mass conserving perturbations in the context of critical dynamics 0, l38j . 

Turning to the correlation function, we proceed as for the EW case and obtain 



-i/ 4 k 4 (t+s-t'-t"; 



(f7(k,f)rK-k,f')>. 



(28) 
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For Gaussian white noise this then yields in 1 + 1 dimensions the expressions 



C(x-y,M) 



D 



E 



(-l)"|x-y| 



2u 



and 



(2^) 2 ^ (2n)!z4 2 " +1)/4 (3-2n) 
L>r(5/4) 



2n+ 1 



(t + a) (3-2nV4_ (t _ fl j(3-an 



)/4 



C(f, S ) 



3?r2^ 



1/4 



(29) 



(30) 



(f + ,s) 3 / 4 - (t - sf' 4 

for the space-time correlation and the autocorrelation functions. Similarly, for spatial correlated noise we have 

DO 

C(x-y,t,«) = ^(-lraW^lx-yl^Kt + s)-^"- 2 ^- 4 )/ 4 - ( t _ s )-(2n-2 P +d-4)/4] (31) 



with an'(/j) 



n=0 

2 2p D r(p)r((l+2«-2p)/4) 



(2n)!(2p-2n+3)7r 1 /2r((l-2p)/2)y 



(l + 2n-2p)/4 



, the autocorrelation being given by the term with n = 0. 



IV. SPACE-TIME SYMMETRIES OF THE NOISELESS EQUATIONS 

We have seen in the previous Section that both models under consideration show dynamical (critical) scal- 
ing behavior. For critical systems it is well know 0, l25j that the multipoint correlators (3(xi, ti, . . . , x n , t n ) := 
(/ii(xi, t\) . . . /i n (x ra , t n )) satisfy a covariance behavior of the kind 



<3(6xi, & z ii 



,bx n ,b z t n ) = b Xl 



XX + ... + X, 



(3(xi , ii . . . , x n , t„) 



with a constant rcscaling factor b and with some numbers Xj. It has been proposed by Henkel [25j to extend this 
dynamical scaling with b constant to rescaling factors which are space and time dependent, i.e. b — ► b(x,t). For 
the special case z = 2 this can be connected to so-called Schrodinger invariance 0> ^3| . In [3 HH the concept of 
Schrodinger invariance has been extended to the case of an arbitrary value of z, but a concrete comparison of the 
space-time quantities predicted by this theory with specific, exactly solvable models with z ^ 2 is still lacking. Here 
we wish to apply this theory - which is called theory of local scale invariance (LSI) - not only to the case z = 2 (i.e. 
the EW case), but also to the MH case with z = 4. 

LSI proposes to derive expressions for the scaling functions of two-time quantities by looking at the space-time 
symmetries of the noiseless equations of motion. In our case the deterministic equations for z = 2 resp. z = 
4 are obtained by dropping the noise term on the right hand side of Eq. resp. @. The case z = 2 then 
yields the free diffusion equation or, equivalently, when going to complex times, the free Schrodinger equation. The 
maximal kinematic group of space-time symmetries which leave the free Schrodinger equation invariant is the so-called 
Schrodinger group |42l l43l l44l |45| , the elements of which transform space and time in the following way: 



Tlx + vt + a 



at + P 
jt + 6' 



*5 - /?7 = 1 



(32) 



where 1Z is a rotation matrix, whereas v, a and a, /?, 7, 6 are real parameters. We write for this also (x', <') = g(x, t) 
and denote the inverse transformation by 9 _1 (x',<'). Under the action of these group elements solutions ^ of the 
Schrodinger equation transform as 



*' = / 3 (. 9 - 1 (x,i))*(. 9 " 1 (x,t)) 



(33) 



where the companion function f g is known explicitly |43j . The generators of the Schrodinger group, which can be 
considered as infinitesimal version of these transformations, form a Lie- Algebra. 

In |25j generators of space-time transformations have been constructed which act as dynamical symmetries on a 
more general deterministic equation, namely on 



-Xd t + 4r<9, 



* = 0. 



(34) 



where z > is a real number. It is easy to see that Equation (|34|l is equivalent to flj for z = 2 resp. to © for z = 4 
when setting V2 = (4A) _1 resp. 1/4 = — (16A) -1 . The generators of the corresponding Lie algebra are explicitely given 
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by (for d = 1) [13 

X_i = -ft (35) 

y_ e = -ft (36) 

X = -tft - irft - - (37) 

z z 

X x = -t 2 d t - -irft - — t - Ar 2 ft 2 " z 
z z 

-2 7l (2 - z)rft x - z - 7l (2 - z)(l - z)d~ z (38) 
= -t0 r - Xzrd 2 r ~ z - 7l z(2-z)ft 1 - 2 (39) 

Operators like ft 1-2 are so-called fractional derivatives, which we recall in the Appendix B. Hereby the quantities x 
and 71 are related by 

x = ^ + ^(2-z). (40) 

As shown in [2^| these space-time symmetries (I35|) - (I39H can be used to fix the form of the two-time response function 
completely. Using all generators and writing r = x — y, one obtains 

Ra(r,t,s) = (t-s)-^ z cj) LJjL, ) («) 

where the index indicates that this is the result for the noise-free theory. The scaling function 4>(u) satisfies the 
fractional differential equation [25| 

(ft + z\udl- z + 2z(2 - z) 7 ift^ 2 ) 4>{u) = 0. (42) 

We have to stress that the scaling function given in [2f| is not the most general solution of this equation. In Appendix 
C we derive this most general solution for any rational z. As shown in the next Section it is this solution which 
permits us to derive the exact expressions for the space-time response and correlation functions in the MH case by 
exploi ting exclusively the space-time symmetries of the deterministic equation J5J. For z = 2 we recover the known 
result (23 

<j)(u) = 0o exp (-Am 2 ) (43) 

where the numerical factor </>o is not fixed by the theory. For the case z = 4 and d = 1 our new result is (see Appendix 
C for the expression obtained in two space dimensions) 

_ / A 4 \ ! ' 4 /3 5 A 4 

- ( x 4 V /2 3 A A . _ /l 3 A 4 \ 

+ ci u 0^2 u +C20-F2 (-,-, u (44) 

1 V 16 J V 4 2' 16 J V 2 4 16 J v ' 

with some constants Co, C\ and C2. Here we have used the fact that x — | in a free field theory as is easily obtained 
from a dimensional analysis (see also the next Section). The coefficients Co, c± and £2 have to be arranged in such a 
way that </>(u) vanishes for u — > 00. Analysing the leadings terms [4rl l47j one realizes that the condition 

r (I) r (!) * + r (!) r (I) * + r (I) r (I) * - » 

provides exactly this, as it cancels all exponentially growing terms. 

Let us close this Section by noting that in the derivation of expression (|41|) we exploited the fact that the exact 
responses of both the EW and MH case are time translation invariant. Often when discussing out-of-equilibrium 
systems this is not the case and one has to consider the sub-algebra where the generator A_ 1; responsible for 
time translation invariance, is omitted jU 0, 13, 113, This then yields for the autoresponse the expression 

Ro(t,s) = ros~ 1 ~ a (j:) 1+a - A «/ 2 (l — l) _1_a , where the parameters a and a' have to be determined by comparing 
with known results. When setting a = a! = Xr/z — 1, one recovers our result. 
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V. DETERMINATION OF RESPONSE AND CORRELATION FUNCTIONS FROM SPACE-TIME 

SYMMETRIES 



In order to use the symmetry considerations of the last Section, we have to adopt the standard field theoretical 
setup for the description of Langevin equations [23l l49l IHcj. Apart from the field h(x, t) wc consider the so-called 
response field h(x, t) which leads to the action 



S[h, h] 



du dR 



h(d u + v 2l {-V 2 ) l )h 



-- / dudRdu' dR'h(R,u)(r](R,u)r](R',u'))h(R',u') 



(46) 



where I = 1 for the EW case and I = 2 for the MH case. The temporal integration is from to oo whereas the spatial 
integration is over the whole space. The following reduction to the noise-free theory is a well known procedure j23| 
but we recall the most important steps in order to establish notations. 
Varying the action yields the equation of motion for the fields h and h 



dh(x, t) 

at 

dh{x, t) 
dt 



-^2;(-V 2 )'/i(x,i) - / dudRh(R,u)(r)(R,u)r)(x,t)) 
v 2l (-V 2 ) l h(x,t) 



(47) 
(48) 



As expected, one recovers for the height h(x, t) the Equations and |J5J by identifying the noise with the term 
- JdudR h(R, u) (7?(R, u)r)(x, £)) . 

One can now proceed by looking at the multi-point functions which are defined in the usual way via functional 
integrals: 



In rn \ r n m 

(Y[h(xi,U) II h(xj,tj)\:= V[h]D[h]Y[h(xi,U) Y[ h(xj,tj)exp(-S[h,h]). 

\i— 1 j—n+1 I i—1 j—n+1 

Within this formalism the space-time response JHJl is given by 

R(*,y,t,s) = (h(x,t)h(y,s)). 
In order to proceed one splits up the action in the same way as done in |23j . that is as 

S[h,h] = So[h,h} + S th [h,h] 

with the deterministic part 



Sn\h,h] 



du dR 



h(R, u) (d u + v 2 i(-V 2 ) l )h(R,u) 



and the noise part 



Sth[h,h] = - I dudRdu 1 dR' h(R,u)(r ] (R,u)r 1 (R , ,u , ))h{R',u') 



(49) 

(50) 
(51) 

(52) 
(53) 



We call the theory exclusively described by So noise-free and denote averages with respect to this theory with (. . .)q. 
The n-point functions of the full theory can then be written as 



Y[h(xi,ti) n ~ h ^^ 

U'=l j=n+l 



'J h(xi,U) h(xj,tj)exp(-S t h[h, h]) 

\i=l j=n+l I 



(54) 



It is easy to see that the noise-free theory has a Gaussian structure both for the EW and the MH model. Introducing 
the two-component field ^ = Cfy one can write the exponential exp(— So[h, h]) as exp(— J du dr du' dr' Aty) with 



2 



6(u - u')6(r - r')((-V 2 )' - d u ) 

6(u-u')S(r-r')((-V 2 ) 1 +d u ) 



(55) 
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From this one deduces two important facts which wc will need in the sequel. Firstly, one has 

(h...hh^Ji)o = (56) 

n m 

unless n = m, which is due to the antidiagonal structure of A (see for instance @, chapter 4) For z = 2 Eq. 1|56|) 
coincides with the Bargmann superselection rule [5lJ. Secondly, Wick's theorem holds. With this it follows that one 
can write the four-point function as 

(h(x,t)h(y,s)h(R,u)h(R,u')) = (h(x,t)h(R,u)) {h(y,s)h(R',u')) (57) 

+ (h(y, s)h(R, u)) (h(x, t)h(R', u')} 

where we have used Eq. (|56[1 . 

Now one can calculate the quantities of interest, namely the space-time response and correlation functions. For this 
one develops the exponential in l|54|) in a power series. One remarks immediately that due to the selection rule 1501) 
one has 

•R(x,y,M) = #o(x,y,i, s) := (h(x,t)h(y,a)) , (58) 

i.e. the linear response function of the full theory is equal to the noise-less linear response function. It follows that 
for any realization of the noise one gets the same expression for the response function, in agreement with the exact 
results derived in Section III. It is also worth noting that within a free field theory non-linear responses vanish due 
to the same superselection rule. Things are of course more tricky for a field theory which is not free as here the noise 
can contribute to the reponse. Whether this is the case depends on the concrete form of the interaction. 

By expanding the exponential in l|54|) we obtain in a similar way that the space-time correlation function is given 
by the expression 

C(x,y,t,s) = J dudKdu' dK' (h(x,t)h(y,s)h(R,u)h(R',u')) (r 1 (Il,u)T](R',u')}. (59) 

Using Wick's theorem we can replace the four-point function by two-point functions (see Eq. 1(57(0 and obtain 

C{x,y,t,s)=2 J dudRdw'dR'(?7(R,w)77(R',u'))(/i(x,t)/i(R,u))o(/i(y,s)/i(R / ,u / ))o. (60) 

Inspection of Eqs. (|58|) and (|60J) reveals that the only remaining undetermined quantity is the two-point function 
(h(x, t)h(y, s))o- However, as duscusscd in the previous Section, this two-point function is fully determined by the 
space-time symmetries of the deterministic equation of motion. It remains to show that the insertion of this two-point 
function into Eqs. 1(58(1 and 160|) indeed yields the exact expressions for the space-time quantities both for the EW 
and for the MH case. 



A. z — 2: The Edwards- Wilkinson case 

This case with Gaussian white noise has already been discussed in p3[ in the context of phase ordering kinetics 
and of critical dynamics. The response function can be read off from the Equation (|41|l after inserting the scaling 
function 1 (43(1 . Recalling that for the EW case we have z = 2, x = 5, and v-x = (4A) _1 we readily obtain the exact 
result 116( 1 . The only quantity left free by the theory is the numerical prefactor </V 

The space-time correlation function is obtained by inserting the same two-point function into Eq. (|60|l . This is 
most easily seen by using the integral representation (|15f) of the two-point function which yields after interchanging 
the order of integration: 

t s 

C(x,y,t,s)=c J dt'J dt"J ^ je ^(x-y) e -, 2 k 2 ( t -O e - 2k 2 (-*") (7}(k;t ' )?K _ kj ^ )) (61) 


which is exactly the same expression (up to the undetermined constant Co) as (|18fl. and this for any choice of the noise 
correlator (f)(k, t')fj(— k, t")). It then immediately follows that we recover the exact results for C given in Section 
IIIA. 
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B. z = 4: The Mullins-Herring case 

For the MH case we proceed along the same lines as for the EW case. Let us start with the one-dimensional case 
d = 1. As already seen for the EW model, the response function is just the response function of the noise-free theory 
whose scaling function is given by l|44|l together with the condition l|45|l needed for a response which vanishes for 
u — > oo. In order to proceed further we remark that in a powers series expansion of (|44|l odd powers of the scaling 
variable u only enter through the term with coefficient Cq. However, odd powers of the scaling variable are absent in 
the exact result l|26[) . so we have to set cq = 0. From Eq. (|45J) it then follows that 




Recalling that for the MH case z = 4, x = |, and v± = — (16(A) _1 , it is now easy to check that the proposed scaling 
function together with (|62[l indeed yields the exact result l|2()|) up to the normalization constant £2 . This so determined 
two-point function can then be inserted into the correlation function l|60|) , yielding the exact result l|28(l . This is again 
most easily seen by using the integral representation 

i?(x - y,t, s) = to J ^- e^-y> e -^(t- s) (63) 

of the response function (|26() . 

In two dimensions we have to replace the expression Q44JI by l|C16(l . It follows that in d = 2 the scaling function 
obtained from LSI only contains two parameters which are furthermore related through the condition I|C17J| . We are 
therefore left with a single undetermined parameter which only appears as a numerical prefactor, similar to the EW 
case. Inserting the resulting scaling function into the expressions (|58|l and (|60|l readily yields the exact results for the 
space-time quantities in two dimensions. 

VI. MICROSCOPIC GROWTH MODELS AND SPACE-TIME CORRELATIONS 

Many theoretical studies of growth processes focus on atomistic models where particles are deposited on a surface 
and are then incorporated into the growing surface following some specific rules which might include local diffusion 
processes. Of special interest is the determination of the universality class to which these models belong. This is 
usually achieved by computing some universal quantities through numerical simulations and comparing them to the 
corresponding quantities obtained from continuum growth equations like the EW and the MH equations discussed in 
this paper. In a commonly used approach one focuses on the estimation of the exponents z and £, which govern the 
behavior of the surface width Q, through the best data collapse. 

In order to show that it is useful to look at two-time quantities in nonequilibrium growth processes we discuss in 
the following the space-time correlation function in the Family model |33j | and in a variant of this model |36j | . Even 
so these are very simple models, there is still some debate on the universality class to which these models belong, 
especially in 2+1 dimensions. Whereas earlier numerical studies yielded the value z = 2 for the dynamical exponent 
in the 2+1-dimensional Family model [^3, 0, E3 , in agreement with the EW universality class with Gaussian 
white noise, Pal et al. [36l l3^ | in their study obtained a value z « 1.65, pointing to a different universality class. In 
addition they studied a variant of this model (which we call restricted Family model in the following) for which they 
recovered z = 2. These results of Pal et al. are surprising, especially so as Vvedensky succeeded j53j in deriving in 
1+1 dimensions the EW equation with Gaussian white noise from both the Family and the restricted Family model 
through a coarse-graining procedure. 

The Family model is a ballistic deposition model with surface diffusion where a particle is dropped at a randomly 
chosen surface site. Instead of fixing itself at this site, the particle first explores the local environment (usually one 
restricts this exploration to the nearest neighbors) and fixes itself at the lattice site with the lowest height. When two 
or more lattice sites other than the originally selected site have the same lowest height, one of these sites is selected 
randomly. In case the originally chosen lattice site is among the sites with the lowest height, the particle remains at 
this site. In the restricted version of this model, introduced in [3(|, the particle only moves to a site of lowest height 
when it is unique. This change has the effect that the moving of the particle only contributes deterministically to the 
surface shape. 

We have simulated these two models both in 1+1 and 2+1 dimensions. For the 1+1 dimensional models all previous 
studies agree that z = 2 and that both models belong to the one-dimensional EW universality class with Gaussian 
white noise. Our main interest here is the height-height space-time correlation function C(x — y, t, s). From the exact 
results presented in the first part of the paper we conclude that this two-time quantity should only depend on the 
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FIG. 1: Dynamical scaling of the space-time correlation function C(r 2 /s, t/s) for Family model in 1+1 dimensions with different 
values of the waiting time s: (a) C vs r 2 /s for some fixed values of t/s, (b) C vs t/s for some fixed values of r 2 /s. The green 
curves are obtained from the exact result I19H derived from the continuum EW equations with uncorrelated Gaussian white 
noise. Numerical error bars are smaller than the sizes of the symbols. 



TABLE I: Estimates for the nonuniversal constants D and V2- 





D 




d = 1 


8.85(4) 


1.260(6) 


d = 1 restricted 


9.25(5) 


1.312(7) 


d= 2 


83(2) 


1.49(5) 


d = 2 restricted 


271(5) 


2.63(6) 



two scaling variables t/s and r 2 / s where r = |x — y|. In Figure 1 we test this expected scaling behavior in the 1+1 
dimensional Family model. In Figure la we fix t/s und plot the correlation function as a function of r /s, whereas 
in Figure lb r 2 / s is fixed and C is plotted vs t/s. Lattices with 12800 sites have been simulated and the data shown 
result from averaging over 1000 runs with different random numbers. The curves obtained for different values of 
the waiting time s collapse on a common master curve when multiplying C with s" 1 / 2 . In addition, these master 
curves nicely agree with the expression l|19(l obtained from the EW equation with uncorrelated white noise, once the 
nonuniversal constants D and have been determined |56j . A similar good agreement is obtained for the restricted 
Family model. We list our estimates for D and v-i for both one-dimensional models in Table [I] It follows from this 
table that D and vi slightly differ in both models. In addition, D/h>2 is slightly larger for the restricted model, even 
so the error bars are overlapping. 

It is worthwhile noting that the continuum description is not expected to completely describe the lattice models for 
small values of r 2 /s, as here the discrete nature of the lattice can not be neglected any more. Deviations are indeed 
observed in Figure lb for r 2 /s ~ 1 and small values of t/s. These lattice effects are expected to get stronger when 
the dimensionality of the system increases. 

After having verified that the computed scaling functions in both versions of the 1+1 dimensional Family model 
agree with the solution of the EW continuum equation, let us now proceed to the more controversial 2+1 dimensional 
case. In Figure 2 we display the space-time correlation computed for the original Family model in 2+1 dimensions. 
Again, in the left panel we fix t/s, whereas in the right panel r 2 /s is kept constant. The data shown here have 
been obtained for lattices with 300 x 300 sites with 5000 runs for every waiting time. We carefully checked that our 
nonequilibrium data are not affected by finite-size effects. Furthermore, we ran different simulations with different 
random number generators and obtained the same results within error bars. We obtain as the main result of these 
simulations that the scaling function of the space-time correlation function is in excellent agreement (once the values 
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FIG. 2: The same as in Figure 1, but now for the Family model in 2+1 dimensions. Numerical error bars are comparable to 
the sizes of the symbols. 



of the nonuniversal constants have been determined, see Table 1) with the exact result obtained from solving the two- 
dimensional EW equation with uncorrelated white noise. The only discrepancies observed for small r 2 / s and small 
t/s are qualitativly the same as for the 1+1 case and are explained by the discrete nature of the lattice. Our results 
arc in accordance with the results of 0, EE EM E3 but strongly disagree with those of Pal et al. |HH E3 ■ Indeed, 
a noninteger value of z in a continuum description can not be realized in a linear stochastic differential equation and 
leads to completely different scaling functions as those obtained from the EW equation. 

In Figure 3 (see also Table 1) wc show our results for the restricted family model in 2+1 dimensions. Again, 
dynamical scaling is observed, and again the data arc well described by the EW scaling functions in the scaling 
limit. However, the determined values of the nonuniversal quantities D and v-i are markedly different from the values 
obtained for the original model. Specifically, the ratio D/v2 (which is of the dimension ksT) is much larger for the 
restricted model. Identifying Djvi with a (nonequilibrium) temperature, we can view the processes in the restricted 
model to take place at a higher temperature than in the original model. This is in agreement with the observation 
from Pal et al. E3 that the surface is locally rougher in the restricted model, as evidenced by the larger value of the 
interface width. In addition, the change in the diffusion rule leads to a nonmonotoncous behavior of the correlation 
function for small r 2 /s, as shown in the inset of Figure 3a for ,s — 25 and t/s = 1.04. Plotting the correlation 
function in both the (10) and the (11) direction, we see that correlations between nearest neighbours are suppressed, 
whereas the autocorrelation, i.e. the correlation with r = 0, is strongly enhanced. This behavior can be understood by 
recalling that in the restricted model a particle only diffuses to a lower nearest neighbor site when this site is unique, 
but otherwise remains on the original site. If we increase s and r, this effect weakens, and it completely vanishes in 
the scaling limit of large waiting times and large values of r 2 / s. 

From our observation that the numerically computed space-time correlation functions of both microscropic models 
coincide in the scaling limit with the exact expression from the EW model we conclude that both the Family model 
and the restricted Family model belong to the EW universality class with uncorrelated noise, and this not only in 
1+1 dimensions but also in 2+1 dimensions. 



VII. CONCLUSIONS 



The aim of the present paper is twofold: on the one hand we discuss the usefulness of universal scaling functions of 
space-time quantities in characterizing the universality class of nonequilibrium growth models, on the other hand we 
demonstrate how a general symmetry principle allows to derive scaling functions of two-point quantities for equilibrium 
and nonequilibrium processes described by linear stochastic Langevin equations. 

In the context of nonequilibrium growth processes it is rather uncommon to study universal scaling functions of two- 
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FIG. 3: The same as in Figure 2, but now for the restricted Family model in 2+1 dimensions. The inset in (a) shows the 
correlation function in the (10) and (11) directions for the case s — 25 and t/s = 1.04. The change of the diffusion rule has a 
strong impact on the autocorrelation with r = and on the nearest neighbor correlations. Numerical error bars are comparable 
to the sizes of the symbols. 



point quantities in the dynamical scaling limit in order to determine the universality class to which a given microscopic 
model belongs. Wc have illustrated the usefulness of this approach by comparing the numerically obtained space-time 
correlation functions for two atomistic growth models with the exact expressions obtained from the corresponding 
continuum stochastic Langevin equation. This approach has allowed us to show that both models belong to the same 
universality class, thus correcting conclusions obtained in earlier numerical studies. 

The study of universal scaling functions of space-time quantities should also be of value in more complex growth 
processes which are no more described by linear stochastic differential equations. Examples include ballistic depo- 
sition with an oblique incident particle beam ^| or growth processes of the KPZ |54| and related j5f| universality 
classes where nonlinear effects can no more be neglected. In addition, the scaling functions studied here can also be 
measured in experiments involving nonequilibrium or equilibrium interface fluctuations. A promising system is given 
by equilibrium step fluctuations [lOL fill Il3. fl9j| . as these are again described by linear Langevin equations. 

In addition we have shown that in nonequilibrium growth processes scaling functions of out-of-equilibrium quantities 
can be derived in a model independent way by exploiting the generalized space-time symmetries of the noiseless part 
of the stochastic equations of motion. We have demonstrated explicitly how to proceed in the case of a rational 
dynamical exponent z, following the general ideas formulated by Henkel a few years ago [2||. In the context of 
nonequilibrium growth processes, the cases z = 2 (EW model) and z = 4 (MH model) are of special interest. The 
case z = 2 has already been studied extensively in the past. For the case 2 = 4, however, we present to our knowledge 
for the first time the derivation of nonequilibrium scaling functions by exploiting the mentioned symmetry principles. 
As these scaling functions are found to agree with the exact expressions derived from the MH equation, wc conclude 
that the postulated space-time symmetries and the proposed way for constructing the scaling functions can also be 
valid for other cases than merely the case z = 2. 

Let us end this paper by a general remark on the applicability of the concept of local scale transformations in the 
context of other out-of-equilibrium processes. The data presented in this paper for the microscopic growth models 
nicely show the limitations of the continuum equations in describing atomistic models. Whereas in the limit of large 
times and large spatial separations the numerically computed correlation functions completely agree with the exact 
results from the continuum equation, notable deviations are observed for small times and small spatial separations. 
These deviations reflect the microscopic details of the models (underlying lattice structure, diffusion rules etc.) which 
are not captured by the continuum model. This sheds an interesting light on an ongoing discussion |48l l57l l58l l59j 
on the applicability of the theory of local scale invariance, which, we recall, permits to derive expressions for scaling 
functions starting from the noiseless part of the continuum equations of motion. Clearly, it has to be expected that the 
so derived scaling functions can not fully describe microscopic models in the short time and short distance limit. As 



14 



it is exactly this limit which is in the centre of the mentioned discussion involving the theory of local scale invariancc, 
it seems advisable to take any observed deviations in this limit between numerical data, obtained from simulations of 
microscopic models, and the theoretically derived scaling functions cum grano salts, as these deviations might only 
reflect the microscopic nature of the models. 



APPENDIX A: EXACT RESULTS FOR THE MULLINS-HERRING CASE IN d = 2 



We here compile the exact results in two space dimensions for the dynamical two-time quantities for the Mullins- 
Hcrring case with z = 4. The space-time response can again be expressed through generalized hypergeometric 
functions and reads 



i?(x-y,t,s) 



For the autoresponse one gets 



8tw, 



1/2 



^2 256i 



l*-y| 



|x-y| 

Au\' 2 \(t-s)^ 



2' 256i/ 4 \(t-sy/ 



-1/2 



(Al) 



(A2) 



As already noted before, these expressions give the response of the system to the noise itself and therefore do not 
depend on the concrete realization of the noise as long as it is nonconserving. 

For the space-time correlation function C(x — y, t, s) and for the autocorrelation function C(t, s) we obtain the 
following expression, which depend on the form of the noise: 

Gaussian white noise (MH1): 



C(x-y,t,s) 



D 
32^ 



y n+1 ' 



~ (-l)»|x-y| 2 »r(? 

(2n)!r(l + 2n)(r(i - n)) 2 *4" +1)/2 (l n) 



C(t,s) = 

Spatially correlated noise (MH2): 



(t + s) 



D 



(l-n)/2 



(l-n)/2 



32^ 5 / 2 ^ 4 1/2 



(t-s) 



(t + s) 1 ' 2 -{t-s) 1 ' 2 



|x-y| 

47T^4 



hi 



t + S 

t-s 



(A3) 
(A4) 



C(x-y,t jS ) = ^(-ira( 2 '(p)|x-y| 2 "[(t + S )-( 2 "- 2 "- 2 )/ 4 
C(M) = 4 2) (p) ^ + S ) (2+2p)/4 -(t- S ) (2+2 " )/4 " 



f2) 

with a„ (p) 



g r( P )r((2+27t-2p)/4) 



2 2 ™(n!) 2 (2p-2n+2)r(l-p)i/ 



(2 + 2n-2p)/4 • 



(i- S ) 



-(2ri-2p-2)/4 



(A5) 
(A6) 



APPENDIX B: ON FRACTIONAL DERIVATIVES 



We list here the most important properties of the fractional derivates as these will be used in Appendix C in the 
derivation of the scaling function of the space-time response. We stress the point that there are several definitions of 
fractional derivatives available, which are not equivalent. However, we need a special type of fractional derivatives. 
We simply quote the most important properties as given in |2fi| and refer the reader to this reference for a more 
thorough introduction. 
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<9" acts on a function f(r) which can be expanded into the form f(r) = X)^e£ fef e + X^^Lo -^n*^™' ( r )- Here E is 
the set E = /iN + A with /j > and A ^ —{^{n + 1) + m + 1), where n, m e N. £W is the n-th derivate of the delta 
function, is then defined by 

i.) d a r {af{r) + (3g(r)) = ad a r f{r) + pd a r g(r) (Bl) 

T(e + 1) 
r(e-a+l)' 



5Ve = FT^-rTT ^"" + ^ ««,e+n+ir(e + 1)<5 ( "> (r) (B2) 



n=0 

-1— n— a 



m.) d r a <5^ (r) = + V S a . m 5^+ m ^ (r) (B3) 

1 —a — n * — ' 

v ; m=0 

Here a and P are real constants and g(r) is another function which can be expanded in the same way as f(r). The 
most important properties of these fractional derivatives are: 

d a r +h f{r) - d?d b J(r) = d b r .d a r f{r) (B4) 

[cVV]/(r) = (d r r-rd r )f(r)=adr 1 f(r) (B5) 

d a r f{ar) = a a d a ar f(ar) (B6) 

d a ar f(r) = a- a d?f(r) (B7) 



APPENDIX C: DERIVATION OF THE SCALING FUNCTION 



In this Appendix we outline the derivation of the scaling function <f)(u) for any rational dynamical exponent z, 
thereby correcting the incomplete result given in |25l |. For notational simplicity we do this in one space dimension. 
For z — 2 and z — 4 we give the results in two space dimensions at the end of this Appendix. 

Our starting point is the equation 1)42(1: 

(d u + aud 2 u - z + bdi-A 0(«) = (CI) 

with a = z\ and b = 2z(2 — 2)71. We write the rational dynamical exponent as z = N + 2, where ./V is the largest 
integer equal or smaller than z. We also assume a^O and give the result for a — at the end, as it can be derived 
in exactly the same way. In a first step we rewrite (ICJ II as 

(d* u + aud u + SW(u) = (C2) 

with *Sf(u) = d^~ z (f>(u). In doing so we have used property (|B4|) of the fractional derivative. It is in fact this step 
which enables us to avoid the negative exponents for the fractional derivatives in equation (|C1|I , which are responsible 
for the incomplete result in |25j. The solution ^f(u) of this equation yields then the scaling function <f>(u) through the 
relation 



<P{u)=dl-^{u). (C3) 

Before doing this, it is instructive to consider (|C2|) for the case z = 4. Indeed, Eq. (|C2|) is then a normal differential 
equation of forth order, so that the solution will be, a priori, a linear combination of four linearly independent 
solutions (i.e. it will contain four free parameters). The method applied in |25j . however, only yields one of these 
linearly independent solutions. 

We solve (|C2|I for u > under the additional assumption that the desired solution is nonsingular for u — > 0. 
Furthermore, we require that the scaling function should drop to zero for u — > 00. We then make the ansatz 

00 

^H = ^c„ U f +s , c ^0 (C4) 

n=0 

and suppose s > — 1 [r30l |. This ansatz is introduced into (|C2(1 and yields, because of cq =/= 0, the recursion relation 



(C5) 
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with £ = p + qN as well as the relation 



where the set £ is given by 



P 

S = h 777*3 m (z E, (C6) 
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f ;= -l,0,...,iV-l, 

\ 0,...,N-1, P = Q 

It then also follows that ci = . . . = Cp+ Q _/v_i = 0. With this we obtain after some algebra: 

oo 

*(«) = ^ c m J] a £™V*+f + m (C8) 



mG£ n— 



with 



(7 



(m) _ 



p/q+m . b/a 



(C9) 



r (nz + 2 + +m + l) r 

The c m are free parameters not fixed by the theory. Finally, the scaling function <fi(u) is obtained from Eq. (|C3|) : 

4>{U) = Cm0 (m) («) (CIO) 



with 



n=0 



where the coefficients fe£ m ^ are given by 



(-az) n T ( £ + 1 + m) T (n + E^±Hi + hi*) 
r ( (n -l)z + f + m + 2J T fs/3±H + b JAJ 

We remark that our final result (|C10() is indeed regular for u — > 0, as 6q™' ) = for m = — 1, . . . , N — 2. This is 
readily seen by recalling that r(Z) = oo for I G —No- As already mentioned in [2^, the radius of convergence is infinite 
for z > 1. 

We also note that the number of free parameters is a priori equal to N if z £ N and TV + 1 else. However, there 
might be cases where some of the independent solutions </>( m ) (it) vanish. 
For completeness let us also quote the result for a = 0. In this case 



T + l + m)T [nz + z + m + l 

&M = (_$)" 1* ^ 2 . (C13) 

r((n+ l)z + l)r((n - 1)2 + | + m + 2J 

The expressions and l|44|l for 0(tt) can be obtaind from ljC10|) by setting 2 = 2 (i.e. p = and TV = 2) or 2 = 4 
(i.e. p = and iV = 4) respectively. For the EW case it is important to note that <^°)(it) vanishes as bn^ = for 
every n which immediately follows from the fact that b = 2z(2 — 2)71 = for 2 = 2. The remaining solution (jy^'{u) 
is then just the exponential function 

(1) (u) = cxp(-Aw 2 ) . (C14) 

For the MH case it is the solution <j)^ which vanishes in the free-field case. Indeed, from Eq. H4()(l we obtain 
= — i by recalling that x = \ . It then follows that the Gamma-function V ( p / q + m + 5Z£ J a l wa ys diverges in the 
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denominator of Eq. (|C12|1 . yielding (j)^{u) ~ for every u. We are therefore left with three independent solutions, 
and after relabelling we obtain the final expression (|44|l : 

/ A A 1/4 /3 5 A 
4>{u) = c -— u vF 2 \ u 



16 J V 4 4 ' 16 



Let us add that the asymptotic behavior of the generalized hypergeometric functions qF 2 is well known |46l l47j. 
Recalling that the scaling function should vanish for u — ► oo, we can exploit this known asymptotic behavior in order 
to derive relations between the parameters c m . For the case z = 4 this then yields the condition (|45|l given in Section 

GYI 

Let us finish this Appendix by quoting the resulting scaling functions <j>(u) in two space dimensions. For the EW 
case with z = 2 we get the same expression (IC14II as for the one-dimensional case. For the MH case with z = 4, our 
calculations yield the expression 

= c ~i f-4« 4 V /2 of 2 (li~^)+£ 2a F 2 (i i,-4« 4 ) (ci6) 

with the additional condition 



16 7 V T 2' 16 7 V 2 16 



c x = -4=&- (C17) 



ACKNOWLEDGMENTS 



We acknowledge the support by the Deutsche Forschungsgcmeinschaft through grant no. PL 323/2. This work 
was supported by the franco-german binational programme PROCOPE. The authors are grateful to Malte Hcnkcl for 
useful comments and helpful discussions. 



P. Meakin, Phys. Rep. 235, 189 (1993). 

T. Halpin-Healy and Y.-C. Zhang, Phys. Rep. 254, 215 (1995). 

A. -L. Barabasi and H. E. Stanley, Fractal Concepts in Surface Growth (Cambridge University Press, Cambridge, 1995). 
J. Krug, Adv. Phys. 46, 139 (1997). 

J. Krug, in Scale Invariance, Interfaces, and Non-Equilibrium Dynamics, edited by A. McKane et al. (Plenum, New York, 
1995), p. 25. 

P. C. Hohenberg and B. I. Halperin, Rev. Mod. Phys. 49, 435 (1977). 

U. C. Tauber, Critical dynamics — A Geld theory approach to equilibrium and non-equilibrium scaling behaviour (Cam- 
bridge University Press), in preparation (see also http: / /www. prrys.vt.edu~tauber/utaeuber.html| |. 
Y.-K. Yu, N.-N. Pang, and T. Halpin-Healy, Phys. Rev. E 50, 5111 (1994). 
S. F. Edwards and D. R. Wilkinson, Proc. R. Soc. London Ser. A 381, 17 (1982). 
M. Giesen, Prog. Surf. Sci. 68, 1 (2001). 

D. B. Dougherty, I. Lyubinetsky, T. L. Einstein, and E. D. Williams, Phys. Rev. B 70, 235422 (2004). 

D. B. Dougherty, C. Tao, O. Bondarchuk, W. G. Cullen, E. D. Williams, M. Constantin, C. Dasgupta, and S. Das Sarma, 

Phys. Rev. E 71, 021602 (2005). 

B. Schmittmann, G. Pruessner, and H. K. Janssen, J. Phys. A 39, 8183 (2006). 

W. W. Mullins, in Metal Surfaces: Structure, Energetics and Kinetics (Am. Soc. Metal, Metals Park, Ohio, 1963), p. 17. 

D. E. Wolf and J. Villain, Europhys. Lett. 13, 389 (1990). 

P.-M. Lam and F. Family, Phys. Rev. A 44, 4854 (1991). 

L. Golubovic and R. Bruinsma, Phys. Rev. Lett. 66, 321 (1991). 

S. Das Sarma and P. Tamborenea, Phys. Rev. Lett. 66, 325 (1991). 

O. Bondarchuk, D. B. Dougherty, M. Degawa, E. D. Williams, M. Constantin, C. Dasgupta, and S. Das Sarma, Phys. Rev. 
B 71, 045426 (2005). 

A. Brri, J. M. Pastor, I. Fernaud, I. Brii, S. Melle, and C. Berenguer, Phys. Rev. Lett. 81, 4008 (1998). 

C. Escudero, Phys. Rev. E 73, 020902(R) (2006). 

F. Family and T. Vicsek, J. Phys. A 18, L75 (1985). 



18 



A. Picone and M. Henkel, Nucl. Phys. B 688, 217 (2004). 

M. Henkel, M. Pleimling, C. Godreche, and J.-M. Luck, Phys. Rev. Lett. 87, 265701 (2001). 
M. Henkel, Nucl. Phys. B641, 405 (2002). 

M. Henkel, A. Picone, and M. Pleimling, Europhys. Lett. 68, 191 (2004). 

M. Henkel and M. Pleimling, Phys. Rev. E 68, 065101 (R) (2003). 

F. Baumann and M. Pleimling, J. Phys. A 39, 1981 (2006). 

F. Baumann, S. Stoimenov, and M. Henkel, J. Phys. A 39, 4095 (2006). 

A. Pelissetto and E. Vicari, Phys. Rep. 368, 549 (2002). 

S. Liibeck, Int. J. Mod. Phys. B 18, 3977 (2004). 

P. Calabrese and A. Gambassi, J. Phys. A 38, R133 (2005). 

F. Family, J. Phys. A 19, L441 (1986). 

P. Meakin and R. Jullien, J. Phys. (Paris) 48, 1651 (1987). 

D. Liu and M. Plischke, Phys. Rev. B 38, 4781 (1988). 
S. Pal and D. P. Landau, Physica A 267, 406 (1999). 

S. Pal, D. P. Landau, and K. Binder, Phys. Rev. B 68, 021601 (2003). 
F. Baumann and M. Henkel, in preparation. 

The response function can also be viewed as the response of the system to the noise itself. 

This assumption is also important in connection with the symmetry based approach presented in this paper. A nonvanishing 
intial average of ft(x, 0) would lead to modifications. 

J. Krug, H. Kallabis, S. N. Majumdar, S. J. Cornell, A. J. Bray, and C. Sire, Phys. Rev. E 56, 2702 (1997). 

C. R. Hagen, Phys. Rev. D 5, 377 (1972). 

U. Niederer, Helv. Phys. Acta 45, 802 (1972). 

S. Lie, Arch. Math. Nat. vid. (Kristiania) 8, 328 (1882). 

H. A. Kastrup, Nucl. Phys. 7, 545 (1968). 

E. M. Wright, J. London Math. Soc. 10 (1935) 287; J. London Math. Soc. 27 (1952), 256, erratum. 
E. M. Wright, Proc. London Math. Soc. 46 (1940) 389. 
M. Henkel, T. Enss, and M. Pleimling, J. Phys. A in pr 



cond- 


mat/0605211 


i, 423 


(1973). 



P. C. Martin, E. D. Siggia, and H. A. Rose, Phys. Rev. 
H. K. Janssen, Z. Phys. B 23, 377 (1976). 
V. Bargmann, Ann. Math. 59, 1 (1954). 
F. D. A. Aarao Reis, Physica A 316, 250 (2002). 

D. D. Vvedensky, Phys. Rev. E 67, 025102(R) (2003). 

M. Kardar, G. Parisi, and Y.-C. Zhang, Phys. Rev. Lett. 56, 889 (1986). 

Z.-W. Lai and S. Das Sarma, Phys. Rev. Lett. 66, 2348 (1991). 

We follow here the procedure used in l27f in the context of phase ordering kinetics. 

M. Pleimling and A. Gambassi, Phys. Rev. B 71, 180401(R) (2005). 

H. Hinrichsen, J. Stat. Mech.: Theor. Exp. L06001 (2006). 

E. Lippiello, F. Corberi, and M. Zannetti, |cond-mat/ 0604321 

One can show that this is justified because of the regularity assumption for u — > 0. It also follows from the regularity 
assumption that a term of the type F n S^ n \u) needs not to be included. 



